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O-l Abstract 

Fisher established an explicit correspondence between the 2-dimensional Ising 
model defined on a graph G and the dimer model defined on a decorated version Q 
of this graph [Fis66]. In this paper we explicitly relate the dimer model associ- 
ated to the critical Ising model and critical cycle rooted spanning forests (CRSFs). 
This relation is established through characteristic polynomials, whose definition 
only depends on the respective fundamental domains, and which encode the com- 
binatorics of the model. We first show a matrix-tree type theorem establishing that 
the dimer characteristic polynomial counts CRSFs of the decorated fundamental 
domain Qi . Our main result consists in explicitly constructing CRSFs of Qi counted 
by the dimer characteristic polynomial, from CRSFs of Gi where edges are assigned 
Kenyon's critical weight function [Ken02]; thus proving a relation on the level of 
configurations between two well known 2-dimensional critical models. 



1 Introduction 

In [Fis66], Fisher established an explicit correspondence between the 2-dimensional 
Ising model defined on a graph G and the dimer model defined on a decorated version Q 
of this graph, known as the Fisher graph of G. Since then, dimer techniques have 
been a powerful tool for tackling the Ising model, see for example the book of [MW73]. 
More recently, in [BdTlOb, BdTlOa], we prove fundamental results for the dimer model 
corresponding to a large class of critical Ising models, by proving explicit formulae for 
the free energy and for the Gibbs measure. 

Critical Ising models we consider are defined on graphs satisfying a geometric property 
called isoradiality. When the underlying isoradial graph G is infinite and Z^-periodic, 
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then so is the Fisher graph Q, and we let Qi = Q jTl? be the fundamental domain. 
The key object involved in the explicit expressions of [BdTlOb] for the free energy 
and the Gibbs measure is the critical dimer characteristic polynomial, whose definition 
only depends on the fundamental domain Qi. This polynomial is a generating function 
for configurations related to super- imposed dimer configurations of Qi, referred to as 
'double-dimer' configurations. By Fisher's correspondence, this implies that the dimer 
characteristic polynomial is a generating function for 'double-Ising' configurations. 

In [BdTlOb], we prove that the dimer characteristic polynomial is equal, up to a con- 
stant, to the critical Laplacian characteristic polynomial of Gi = G/1?, where edges 
of Gi are assigned Kenyon's critical weight function [Ken02]. Using a generalization 
of Kirchhoff's matrix tree theorem due to Forman [For93], the latter is shown to be a 
generating function for cycle rooted spanning forests (CRSFs) of Gi. This suggests the 
existence of an explicit relation between 'double-Ising' configurations and CRSFs, which 
we were not able to find in [BdTlOb]. The first result of this paper is a matrix-tree 
type theorem, proving that the critical dimer characteristic polynomial is a generating 
function for CRSFs of the Fisher graph Qi, see Theorem 7 of Section 4.2. Then, the 
main result of this paper can loosely be stated as follows, refer to Theorem 17 of Section 
4.4 and Theorem 29 of Section 6.1 for more precise statements. 

Theorem 1. Consider a critical Ising model defined on an infinite, 1? -periodic isora- 
dial graph G. Then, there exists an explicit way of constructing CRSFs of Qi counted by 
the critical dimer characteristic polynomial, from CRSFs of Gi counted by the critical 
Laplacian characteristic polynomial. 

This exhibits an explicit relation, on the level of configurations, between two well known 
models of statistical mechanics: the Ising model and CRSFs at criticality. Note that 
such a relation was already suspected by Messikh [Mes06]. Before giving an outline of 
the paper, let us make a few comments. 

• The main contribution of this paper is the proof of Theorem 1, where we actually 
provide the explicit construction. 

• The partition functions (weighted sum of configurations) of both models can be 
expressed from their respective characteristic polynomials, so that it is actually 
stronger to work with characteristic polynomials. 

• Working with graphs embedded on the torus has the advantage of avoiding bound- 
ary issues, but has the additional difficulty of involving the geometry of the torus, 
with non-trivial cycles occurring in configurations. Working on finite pieces of 
infinite graphs, and precisely specifying boundary conditions, would certainly ex- 
plicitly relate double-dimer configurations and spanning trees. 

• Spanning trees are a well suited object for defining a height function and prove 
Gaussian fiuctuations. Thus, it might be that Theorem 1 could be used to prove 
results which are numerically described in the paper [Willi]. 
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Outline of the paper 



Section 2: Definition of the critical Ising model and of the dimer model. Description of 
Fisher's correspondence relating the two. 

Section 3: Definition of the critical dimer and Laplacian characteristic polynomials. Relation 
between the Laplacian characteristic polynomial and CRSFs. 

Section 4: Statement and proof of Theorem 7 establishing that the critical dimer charac- 
teristic polynomial is a generating function for CRSFs of the Fisher graph Qi. 
Precise statement of Theorem 1. 

Section 5: Definition and properties of licit primal/dual edge moves, which are one of the 
key ingredients of the correspondence. 

Section 6: Explicit construction of CRSFs of Qi from CRSFs of Gi and proof of Theorem 1. 



2 Critical Ising model and dimer model 

In this section, we define the 2-dimcnsional critical Ising model, the dimer model and 
describe Fisher's correspondence relating the two. 

2.1 Critical Ising model 

Let G = {V(G), E{G)) be a finite, planar, unoriented graph, together with a collection 
of positive real numbers J = {Je)eeE(G) indexed by the edges of G. The Ising model 
on G with coupling constants J is defined as follows. A spin configuration cr of G is a 
function of the vertices of G with values in {—1, +1}. The probability of occurrence of 
a spin configuration a is given by the Ising Boltzmann measure, denoted Pising^ 

IPlsing(o-) = y^exp j ^ Je(7u<7v 
^^^^^g \e=uveE{G) 

where Zigmg = Saei-i.ij^CG) exp {^e=uveE{G) Je<^u(^v^ , is the Ising partition function. 

We consider Ising models defined on a class of embedded graphs which have an addi- 
tional property called isoradiality. A graph G is said to be isoradial [Ken02], if it has 
an embedding in the plane such that every face is inscribed in a circle of radius 1. We 
ask moreover that all circumcenters of the faces are in the closure of the faces. From 
now on, when we speak of the graph G, we mean the graph together with a particular 
isoradial embedding in the plane. Examples of isoradial graphs are the square and the 
honeycomb lattice. Refer to Figure 1 (left) for a more general example of isoradial 
graph. 
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To such a graph is naturally associated the diamond graph, denoted by G^: vertices 
of consist in the vertices of G and the circumcenters of the faces of G (which are 
also the dual vertices of G*); the circumcenter of each face is then joined to all vertices 
which are on the boundary of this face, see Figure 1 (center). Since G is isoradial, all 
faces of are side-length- 1 rhombi. Moreover, each edge e of G is the diagonal of 
exactly one rhombus of G^; we let 9e be the half-angle of the rhombus at the vertex it 
has in common with e, see Figure 1 (right). 



Figure 1: Left: example of isoradial graph. Center: corresponding diamond graph. 
Right: rhombus half-angle associated to an edge e of the graph. 

The same construction can be done for infinite and toroidal isoradial graphs, in which 
case the embedding is in the plane or in the torus. 

It is then natural to choose the coupling constants J of the Ising model to depend on 
the geometry of the embedded graph: let us assume that Jg is a function of 6e, the 
rhombus half-angle assigned to the edge e. 

We impose two more conditions on the coupling constants. First, we ask that the 
Ising model on G with coupling constants J as above is Z -invariant, that is, invariant 
under star-triangle transformations of the underlying graph. Next, we impose that the 
Ising model satisfies a generalized form of self-duality. These conditions completely 
determine the coupling constants J, known as critical coupling constants: for every 
edge e of G, 



The Z-invariant Ising model on an isoradial graph with this particular choice of coupling 
constants is referred to as critical Ising model. This model was introduced by Baxter 
in [Bax86]. A more detailed definition is given in [BdTlOb]. 

2.2 Dimer model 

Let Q = {V{Q),E{Q)) be a finite, planar, unoriented graph, and suppose that edges of 
Q are assigned a positive weight function v = (fe)eG-E(e)- The dimer model on Q with 
weight function v is defined as follows. 
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A dimer configuration M ofQ, also called perfect matching, is a subset of edges of Q such 
that every vertex is incident to exactly one edge of M. Let M.{Q) be the set of dimer 
configurations of the graph Q. The probability of occurrence of a dimer configuration 
M is given by the dimer Boltzmann measure, denoted Pdimer^ 

P IM\ — HegM 
iTdimeri^Wj — — , 

^dimer 

where ^dimer = ^MeM{Q) riegM is the dimer partition function. 



2.3 Fisher's correspondence 

Fisher's correspondence [Fis66] holds for a general Ising model defined on a finite graph 
G embedded on a surface without boundary, with coupling constants J. We use the 
following slight variation of the correspondence. 

The decorated graph, on which the dimer configurations live, is constructed from G 
as follows. Every vertex of degree A; of G is replaced by a decoration consisting of 3fc 
vertices: a triangle is attached to every edge incident to this vertex, and these triangles 
are joined by edges in a circular way, see Figure 2 below. This new graph, denoted 
by Q, is also embedded on the surface without boundary and has vertices of degree 3. 
It is referred to as the Fisher graph of G. 




Figure 2: Left: a vertex of G with its incoming edges. Right: corresponding decoration 
of 

Fisher's correspondence uses the high temperature expansion of the Ising partition 
function, see for example [Bax89]: 

Rising = I n cosh(Je) I 2l^(G)l ^ []tanh(Je), 

\egB(G) / CeV eeC 

where V is the family of all polygonal contours drawn on G, for which every edge of 

G is used at most once. This expansion defines a measure on the set of polygonal 
contours V of G: the probability of occurrence of a polygonal contour C is proportional 
to the product of the weights of the edges it contains, where the weight of an edge e is 
tanh(Je). 
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Here comes the correspondence: to any contour configuration C coming from the high- 
temperature expansion of the Ising model on G, we associate 2^^^^^^ dimer configu- 
rations on Q: edges present (resp. absent) in C are absent (resp. present) in the 
corresponding dimer configuration of Q. Once the state of these edges is fixed, there 
is, for every decorated vertex, exactly two ways to complete the configuration into a 
dimer configuration. Figure 3 below gives an example in the case where G is the square 
lattice Z^. 



1 2 



3 4 




Figure 3: Polygonal contour of 7?, and corresponding dimer configurations of the 
associated Fisher graph. 



Let us assign, to an edge e of G, weight Vg = 1; if it belongs to a decoration; and weight 
fe = coth Je, if it corresponds to an edge of G. Then the correspondence is measure- 
preserving: every contour configuration C has the same number (2l^(<^)l) of images by 
this correspondence, and the product of the weights of the edges in C, f^^g^ tanh( Jg) 
is proportional to the weight He^C °f corresponding dimer configu- 

rations for a proportionality factor, neGE(G) tanh(Je), which is independent of C. 

As a consequence of Fisher's correspondence, we have the following relation between 
the Ising and dimer partition functions: 

■Rising = I JJ sinh(Je) j Zdimer- 

Fisher's correspondence between Ising contour configurations and dimer configurations 
naturally extends to the case where G is an infinite planar graph. 

2.4 Critical dimer model on Fisher graphs 

Consider a critical Ising model defined on an isoradial graph G embedded in the torus, 
or in the plane. Then, the dimer weights of the corresponding dimer model on the 



6 



Fisher graph Q are: 




if e belongs to a decoration, 
(^e) = cot (^) if e comes from an edge of G. 



(1) 



We refer to these weights as critical dimer weights, and to the corresponding dimer 
model as critical dimer model on the Fisher graph Q. 

3 Critical characteristic polynomials 

In this section we define the critical dimer and Laplacian characteristic polynomials. 
We then state Forman's theorem proving that the Laplacian characteristic polynomial 
is a generating function for CRSFs of the underlying graph. 

3.1 Critical dimer chciracteristic polynomial 

The dimer model has the specific feature of having an explicit formula for the partition 
function due to Kasteleyn [Kas61] and independently to Temperley and Fisher [TF61]. 
It involves a weighted adjacency matrix of the underlying graph known as a Kasteleyn 
matrix. Let us define it for the critical dimer model on a Fisher graph Q, which we 
assume to be the Fisher graph of an infinite, Z^-periodic isoradial graph G. Recall that 
Q\ = QI1? denotes the fundamental domain of Q. 

A Kasteleyn orientation of Q is an orientation of the edges of Q such that all elementary 
cycles are clockwise odd, i.e. when traveling clockwise around the edges of any elemen- 
tary cycle of G, the number of co-oriented edges is odd. When the graph is planar, such 
an orientation always exists [Kas67]. For later purposes, we need to keep track of the 
orientation of the edges of Q. We thus choose a specific Kasteleyn orientation of Q in 
which every triangle of every decoration is oriented clockwise. Having a Kasteleyn ori- 
entation of the graph Q then amounts to finding a Kasteleyn orientation of the planar 
graph obtained from Q by contracting each triangle to a single vertex, which exists by 
Kasteleyn's theorem [Kas67]. Refer to Figure 4 for an example of such an orientation 
in the case where G = 1?. 



The Kasteleyn matrix corresponding to such an orientation is an infinite matrix, whose 
rows and columns are indexed by vertices of Q, defined by: 



K, 



= e. 
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Figure 4: An example of Kasteleyn orientation of the Fisher graph of Z^, in which 
every triangle of every decoration is oriented clockwise. 



where 

1 if X ~ y, and x 



^x,y — \ 



— 1 if X ~ y, and x <— 
else, 



and u is the critical dimer model weight function of Equation (1). Note that K can be 
interpreted as an operator acting on C^^^^: 

V/GC^(^), {Kf),= ^ K,,yfy. 

xev{g) 

The critical dimer characteristic polynomial, denoted by -Pdimcr(z, w), is the determinant 
of the Fourier transform of the Kasteleyn operator. More explicitly, let 7h and 7v be two 
paths in the dual graph of Qi winding once around the torus horizontally and vertically 
respectively. Then, the Fourier transform of K is the modified weight Kasteleyn matrix 
Ki{z,\n) whose lines and columns are indexed by vertices of Qi, and whose coefficients 
are those of K multiplied by z^^ (resp. w^^) when the corresponding edge is crossed by 
the horizontal cycle (resp. vertical cycle), and the sign ±1 is defined by the Kasteleyn 
orientation of the edge. The critical dimer characteristic polynomial then is: 

^'dimer(z,w) = detKi(z,w). 

It is the key ingredient used in explicit formulae for the critical dimer model defined on 
the infinite graph Q or on the finite toroidal graph Qn = G/nZ^, n£W, sec [BdTlOb]. 
More precisely, the partition function of the fundamental domain Qi can be expressed as 
a linear combination of the square root of -Pdimer(z, w) evaluated atz, w € {±1}, and the 
partition function of Qn can be expressed using Pdimer(z, w) evaluated at z", w" G {±1}- 
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By expanding Pdimer(z, w), see also the proof of Lemma 13 of [BdTlOb], 

fdimer(z,w)= J] Sgn((7) JJ (iiTi (z, w))^^^(^) , 

(^e«S|v(gi)| xeViGi) 

and using the fact that Ki{z,w) is an adjacency matrix, one observes that the only 
contribution to the sum comes from configurations which arc unions of disjoint cycles 
covering all vertices of Qi, such that trivial cycles (homotopic to a point) are of even 
length (when the length is 2, it is then a doubled edge), non-trivial cycles (with non- 
trivial homology) can be of even or odd length, and each non-trivial cycle contributes a 
term z^w^ where (h, v) is its homology class. Moreover, since cycles are disjoint, non- 
trivial cycles must be parallel. The difference between these configurations and super- 
imposition of dimer configurations, also known as double- dimers, lies in the terms z, w 
and in the fact that non-trivial cycles can be of odd length; double-dimer configurations 
can be recovered by taking a linear combination of Pdijncr(z,w) with z,w E {il}- We 
refer to configurations counted by -Pdimer(z, w) as 'double-dimer' configurations. 

3.2 Critical Laplacian characteristic polynomial 

A generalization of Kirchhoff 's matrix tree theorem due to Forman [For93] proves that 
the Laplacian characteristic polynomial is a generating function for cycle rooted span- 
ning forests, which are the natural pendent of spanning trees when working on the 
torus. In this section we first define cycle rooted spanning forests, then the Laplacian 
characteristic polynomial, and finally state Forman's theorem. 

We let G be an infinite, Z^-periodic isoradial graph, and Gi = Gjl? be the fundamental 
domain. Note that the content of Sections 3.2.1 and 3.2.2 holds in more generality, i.e. 
when G\ is any graph embedded on the torus and p (see below) is any positive weight 
function on unoriented edges of G\ . 

3.2.1 Cycle rooted spanning forests 

A cycle-rooted tree (CRT) of a toroidal graph Gi is a connected subgraph of Gi with 
a unique non-trivial cycle. A cycle-rooted spanning forest (CRSF) is a collection of 
disjoint cycle-rooted trees covering every vertex of Gi, thus implying that all non- 
trivial cycles are parallel. An oriented CRT (OCRT) is a CRT in which edges of the 
branches are oriented towards the non-trivial cycle, and the non-trivial cycle is oriented 
in one of the two possible ways. An oriented CRSF (OCRSF) is a CRSF consisting of 
OCRTs. 

Let us denote by T a generic OCRT of Gi, by F a generic oriented OCRSF of Gi, and 
by T{Gi) the collection of OCRSFs of Gi. 

Remark 2. 
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• To a CRSF naturally corresponds 2l"°°-*"^'^' ^y^'^'^l OCRSFs. 

• A CRSF is characterized as a subset of |F(Gi)| edges of Gi containing no trivial 
cycle. 

• An OCRSF is characterized as a subset of oriented edges of Gi such that each 
vertex has exactly one outgoing edge of this subset, and which contains no trivial 
cycle. 

Let 7h and 7v be two paths in the dual graph of Gi winding once around the torus 

horizontally and vertically respectively. Assume that 7h and 7v are assigned a reference 
orientation. The homology class of an OCRT T, denoted by H{T) = {h{T),v{T)), is 
defined to be the homology class of its non-trivial cycle in 1?. Define the reference 
number of T to be: 

Ho{T) = {ho{T),vo{T)) = ±{h{T),v{T)), 

where the sign is chosen so that, ho > 0, and vq > when Iiq = 0. Note that this 
definition is independent of the orientation of the non-trivial cycle, so that it also 
makes sense for CRTs. Define the sign of the non-trivial cycle of T to be: 

^{T) = l^H{T)=HoiT)} - l{i?(r)=-i?o{T)}- 
Then, the homology class of the OCRT T can be rewritten as: 

H{T) = N{T)Ho{T). 

Let F be an OCRSF of Gi, and denote by Ti, ■ ■ ■ , r„, its tree components, then the 
homology class H{F) of F, is naturally defined by: 

H{F) = {hiF),viF)) = Y.H{T,) = {j2h{Ti),Y.<Ti)). 

1=1 i=l i=l 

Since non-trivial cycles of the CRT components of F are parallel, we deduce that the 
number HQ{Ti), i G {!,••• ,n}, is independent of i. It is then natural to define the 
reference number of the OCRSF F by: 

Ho{F) = {ho{F),vo{F)) = Ho{Ti). 

As a consequence, the homology class of the OCRSF F can be rewritten as: 

H{F) = N{F)Ho{F), 

where N{F) = X;"=i N{Ti) is the signed number of cycles of the OCRSF F. 
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3.2.2 Critical Laplacian characteristic polynomial 



Suppose that (unoriented) edges of G are assigned Kenyon's critical weight function for 
the Laplacian [Ken02], denoted by p, 



where 9e is the rhombus half-angle of the edge e. Then, the critical Laplacian A on G, 
is represented by the following matrix, also denoted A, whose lines and columns are 
indexed by vertices of G: 



The critical Laplacian characteristic polynomial, PLap(z,w), is the determinant of the 
Fourier transform of the Laplacian operator, that is, i^ap(z,w) = det Ai(z,w), where 
Ai (z, w) is the modified weight Laplacian matrix defined in a way similar to the modified 

weight Kasteleyn matrix. 

A remarkable fact due to Kirchhoff is that, when the graph is finite and embedded in 
the plane, the absolute value of any cofactor of the Laplacian matrix yields the weighted 
number of spanning trees. Forman generalized this result, and for the case of the torus, 
his result can be stated as, see also Lemma 9 of [BdTlOb]: 

Theorem 3 ([For93]). The critical Laplacian characteristic polynomial is the following 
combinatorial sum: 



Note that the weight function p is independent of the orientation of the edges. 

4 Matrix-tree theorem for the Kasteleyn matrix 

Consider a critical Ising model defined on an infinite, Z^-periodic isoradial graph G, 
and let Q be the Fisher graph of the corresponding critical dimer model. In Section 4.2, 
we state and prove Theorem 7, which is a matrix-tree type theorem for the Kaste- 
leyn matrix, thus establishing that the critical dimer characteristic polynomial can be 
rewritten as a generating function for OCRSFs of Qi. In Section 4.3, we analyze this 
polynomial and show that the contribution of some OCRSFs cancel out, leading to the 
definition of essential OCRSFs of Qi. Then, in Section 4.4, we give a precise statement 
of Theorem 1, and an idea of the proof. Section 4.1 is dedicated to notations. 



Ve G E{G), pe = p{9e) = tan 6", 



A. 



^x,y — 



' Pxy if X ~ y 

< - Ej/~a; Pxy if a; = 2/ 
otherwise. 
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4.1 Notations 



We specify notations and terminology for the isoradial graph G and for the correspond- 
ing Fisher graph Q, which were introduced in [BdTlOa], see also Figure 5. These will 
be used throughout the remainder of the paper. 

Edges of Q corresponding to edges of G are referred to as long edges, and edges of the 
decorations of Q are referred to as short ones. 

Vertices of the graph G are written in boldface, and those of G with normal symbols. 
Edges and edge subsets of G are also written in boldface. Let x be a vertex of Q, then 
X belongs to the decoration corresponding to a unique vertex x of G. We shall also 
denote by x the decoration in Q. Conversely, vertices of a decoration x of ^ are labeled 
as follows. Let dx be the degree of the vertex x in G, then the decoration xoi Q consists 
of (ix triangles, labeled ti(x), • • • , t,i^(x) in counterclockwise order, and dx inner edges. 
Vertices of the triangle tk{x) are labeled t'fc(x), Wfc(x), 2;fc(x) in counterclockwise order, 
where t'fc(x) is the only vertex incident to a long edge. We also refer to the triangle 
ffe(x) as the triangle of the vertex Vk{yi)- Whenever no confusion occurs, we drop the 
argument x in the notations above. 



Figure 5: Left: vertex x of the graph G. Right: corresponding decoration x of Q. 

Define a vertex x of to be of type 'v if x = (x) for some vertex x of G and some 
k & {!,■ ■ ■ , dx}, and similarly for 'it;' and 'z'. 

The isoradial embedding of the graph G fixes an embedding of the corresponding di- 
amond graph G^. There is a natural way of assigning rhombus unit- vectors of to 
vertices of Q: for every vertex x of G, and every k G {I,-- - ,dx}, let us associate 
the rhombus unit- vector e^'^^k^^) to Wk{x), e*°^fc(''^ to Zk{x.), and the two rhombus- unit 
vectors e*""'*^^), e*"^fc(^) to Vfe(x), as in Figure 6 below. Note that e*""'fe(^) = e*""^'=+i'='^ 

4.2 Matrix-tree theorem for the Kasteleyn matrix 

In this section, we prove a matrix-tree type theorem for the Kasteleyn matrix Ki{z,w) 
of the critical dimer model on the graph Qi. A key requirement for such a theorem to 
hold is to have a vector in the kernel of the matrix K, which is the subject of the next 
section. 





inner edge 



long edge 
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G and G 




Figure 6: Rhombus vectors of the diamond graph assigned to vertices of Q. 



4.2.1 Vector in the kernel of the Kasteleyn matrix 

The vector in the kernel of the Kasteleyn matrix K is naturally obtained by setting 
A = in the definition of the complex-valued function introduced in Section 4.2.2 of 
[BdTlOa]. More precisely, we define / = (/x){xeV(e)}i by: 



fx 



. "tt.fc(x) 

e * 2 if X = Wfc(x) 

— e~* 2 \i X = 2fc(x) (2) 

, fwk (x) + fzk (x) if X = Wfe (x) , 



for every vertex x of G, and every k G {1, • • • , dx}. Then, setting A = in Proposi- 
tion 15 of [BdTlOa] yields: 

Lemma 4. The vector f is in the kernel of the matrix K. That is, if we let x be a 
vertex ofQ, and xi,X2,x^ be its three neighbors, then: 



{Kf), = Y,K.,xJx, = ^- 



1=1 

Remark 5. In order for the vector / to be well defined, the angles att,j.(x)i Oiz^.{x) need 
to be well defined mod 47r, indeed half-angles need to be well defined mod 2tt. The 
latter are defined inductively in [BdTlOa] as follows, see also Figure 7. Note that the 
definition relies on our choice of Kasteleyn orientation of Section 3.1. Fix a vertex xq 
of G, and set a2;^(xQ) = 0. Then, for vertices of G in the decoration of a vertex x € G, 
define: 

Q?«)j.(x) = '^2j.(x) + 2^A;(x), where 0fc(x) > is the rhombus half-angle of Figure 7, 

aw^{x) if the edge Wfc(x)zfc+i(x) is oriented from Wki'x.) to Zfc_|_i(x) 

a«;fc(x) + 27r else. 

(3) 



a 



2fc + l(x) 



13 



Here is the rule defining angles in the neighboring decoration, corresponding to a vertex 
y of G. Let k and i be indices such that i'fc(x) is adjacent to ^^(y) in Q. Then, define: 




Figure 7: Notations for the definition of the angles in R/47rZ. 

Lemma 6 ([BdTlOa]). For every vertex li. of G , and every k {1, ■■ ■ ,d(x)}, the angles 
^Wf,(x)! CK2fe(x); ^'^s '^^^^ defined in M/47rZ. 



4.2.2 Matrix-tree theorem 

The matrix-tree type theorem is most easily written for the matrix K^, which is the 
following gauge transformation of the Kasteleyn matrix K: 

= D*KD, 

where D is the diagonal matrix whose elements are indexed by vertices of Q, and such 
that Dx,x = fx- Let P^imeii^^^) ~ det Ki{z,\n) be the characteristic polynomial of the 
matrix K^, then we clearly have: 

^'dimer(z,w)= ( J] 1/^ N ^dimer (z, w) . (5) 

Theorem 7. [Matrix-tree theorem for the Kasteleyn matrix] 

The critical dimer characteristic polynomial ^'dimei.(z, w) is the following combinatorial 
sum: 

^d°n.er(z,w)= Y: i n fxTyKxAYlil-Z^^^-^^^^^^ 

FGTiGi) \e={x,y)€F J TeF 
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In order to prove Theorem 7, we introduce two modified incidence matrices M and N 
associated to the graph Q, defined as follows. Rows of M are indexed by vertices of Q, 
columns by unoriented edges of Q, and: 



Kx^y if the edge e is incident to x, and e = xy 
if e is not incident to x. 



Rows of N are indexed by unoriented edges of Q, columns by vertices of Q, and: 



N 



fxfy if the edge e is incident to x, and e = xy 
if e is not incident to x. 



The next lemma relates the Kasteleyn matrix K^, and the incidence matrices M and N. 
Lemma 8. The following identity holds: 

= MN. (6) 

Proof. Let x and y be two vertices of Q. If x and y are at distance more than 2, then 
clearly (MN)x,y = = K^^y. If x and y are neighbors, we have: 

(MAT),,, = Mx,eNe,y = Kx,yfyTx = {D*KD),,y = Kly. 

li x = y, then: 

(MiV)^,^ = J2 M^,eNe,x = Yl ^x,yfxh 
e incident to x 

= 0, (by Lemma 4) 



x.x 



□ 



Proof of Theorem 7. The proof is similar to that of the matrix-tree theorem on the 
torus, see [BdTlOb] for example. Taking the Fourier transform of Equation (6) yields, 
K^{z,w) = Mi(z2, w2)A''i(z2,w2), where Mi(z,w) and Ni{z,\n) are the Fourier trans- 
form of the matrices M and N. Moreover, since for (z,w) € the matrix K^{z,v\i) is 
skew-hermitian, we have: 

^dimer(z,w) = dct (z-^ , w"^ )iVi (z-^ Vy-^ )) , 
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where we also use the fact that = 6|£'(Gi)| is even. We now use Cauchy-Binet's 

formula to expand ^^{^^^^.(z, w): 

^'dimer(z,w) = ^ det(Mi (z'^ , ^ )S) det (iVi (z'^ , w" 5 . (7) 

s c E{g,) 
\s\ = \vig,)\ 

Recall that 7h, 7v are two paths in the dual graph of winding once around the torus 
horizontally and vertically respectively. Suppose first that S contains a trivial cycle, 
denoted by C, which does not cross the horizontal or the vertical cycle 7h, 7v Let us 
show that in this case, det(Mi(z~2, w~2)'S') = o. Define (p : V{Qi) -)■ C by: 



fx 



1 if X G C 
else. 



Then, for every unoriented edge e G 5, we have: 

_Jo ife^C 

\fxMi{z~^,\N-^)x^e + 'PyMi{z-^,\N-^)y^e = K^,y + Ky^^ = if e = G C. 



Thus, det(Mi(z-2,w-^)^) = 0. If C crosses the horizontal and/or vertical cycles 7h, 
7v, then the vector ip in the kernel of (Mi(z^2 ^ w^a)'^')* can be defined in a similar 
way. Therefore, the only contribution to the sum (7) comes from graphs whose number 
of edges equals the number of vertices and which contain no trivial cycle, i.e. from 
CRSFs. 



Let us now compute the contribution of a CRSF F. After a possible reordering of 

11 11 
the rows and columns of Mi(z^2jW^2) and A'^i(z^2 , w^a we can suppose that both 

these matrices are block diagonal, each diagonal block corresponding to a connected 

component of F, i.e. a cycle rooted tree. 

The determinant of a CRT T in Mi(z~2,w^2) (resp. A^i(z^2 , w^2)) can be evaluated 
by expanding it along columns (resp. rows) corresponding to leaves of the CRT. What 
remains then is the evaluation of the determinant reduced to the cycle. More precisely, 
suppose that edges of the branches are oriented from the leaves to the non-trivial cycle. 
Then, the contribution of the branches to -Pdimer(^' 

Kx,yfxfy 

{x,y)E branch of T 

Recall the definition of the reference number Ho{T) = {ho{T),vo{T)) of the CRT T, 
given in Section 3.2.1, and let xi, • • • , x„, be a labeling of the vertices of its non-trivial 
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cycle in the direction given by Hq(T). Then, the contribution of T to Mi[z 2 ^ w 2 ) is: 

I lli^a.„Xi+i 1 Z 2 w 2 +(-lf+M 1 z 2 w 2 , 

and the contribution to Ni[z~2 ^w~^) is: 

I 1 Z 2 w 2 I 1 z 2 w 2 . 

Since the Kasteleyn matrix K is skew-symmetric, we have (— l)""*"^ niLi -^a;<+i,a:i ~ 
~ nr=i -^3;i,a;i+n thus the contribution of the CRT T is: 

{p^Uj^.K^.,..,)l (l-z'^o(T)^.o(T))+ (l-z-'^o(T)^-.o(T)). 

Let Ti and r2 be the two OCRTs corresponding to T, then the contribution of the 
CRT T can be rewritten as: 

j=l \e={x,y)eTi j 

Taking the product over ah CRTs of the CRSF F, and summing over all CRSFs yields: 

pi^.A-M= E n e( n /.^O (i-z'^(^^v(^')), 

Fe{CRSF of Qx) T6{CRT of F} i=l \e={x,y)^Ti j 

E nf n «0(i-z^(^v(^)). 

FeT{gi)TeF \e={x,y)eT J 



□ 



4.3 Characterization of OCRSFs contributing to i^dimer(^5 ^) 

In this section, we characterize OCRSFs of Gi which contribute to -Primer Indeed, 
it turns out that the contributions of some of them cancel out. 

More precisely, let L be a subset of oriented edges of Gi, defining a subset of oriented 
long edges L of ^1. An OCRSF F of Qi is said to be compatible with L, if the long 
edges of F are exactly those of L. We first characterize OCRSFs compatible with L, 
and then OCRSFs compatible with L which actually contribute to P °j^ p^(z, w). 
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4.3.1 OCRSFs compatible with L 



An oriented edge (x, y) of L corresponds to a unique oriented edge (^^(x), ^^(y)) of Qi. 
We refer to tlie vertex ffc(x) as a root vertex of L. This defines for every decoration 
X, a set of root vertices denoted by i?x(L), and a set of non-root vertices of type 'i;', 
i?x(L)'^ := {^i(x), • • • , \ i?x(L), see Figure 8 below for an example. 




Figure 8: In this example L is drawn in bold lines, i?x(L) = {v2(x), ti3(x)} and 
R^{Ly = {vi{^),V4{^)}. 

Remark 9. Suppose that one of the decoration x of Gi contains no root vertex of L, 
and suppose that there is an OCRSF F of Gi compatible with L. Then, since by 
Remark 2, the OCRSF F has exactly one outgoing edge at every vertex, this implies 
that F restricted to the decoration x has a number of edges equal to the number of 
vertices of the decoration, thus it must contain a trivial cycle and cannot be an OCRSF. 
As a consequence, we only consider subset of oriented edges L of Gi such that: 

Every decoration of Gi has at least one root vertex of L. (*) 

Lemma 10. A subset of oriented edges F of Gi is an OCRSF compatible with L, iff 
the following conditions hold: 

1. long edges of F are exactly those ofL, 

2. F contains no cycle consisting of long and short edges, 

3. the restriction of F to every decoration x, is an oriented spanning forest with 
|i?x(L)| connected components whose roots are distinct elements o/i?x(L). 

Remark 11. Condition 3. is equivalent to Condition 3': the restriction of F to every 
decoration x contains no trivial cycle, and is such that vertices of i?x(L) have no 
outgoing edge of F, and every other vertex has exactly one outgoing edge of F. 

Proof. Let F be a subset of oriented edges of Gi- Then by the geometry of the graph 
^1, F cannot have trivial cycles consisting of long edges only. Thus, Conditions 1, 2, 3' 
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are equivalent to saying that the oriented edge configuration F is compatible with L, 
contains no trivial cycle, and is such that every vertex of Qi has exactly one outgoing 
edge of F. By Remark 2, this is equivalent to saying that F is an OCRSF compatible 
with L. □ 

4.3.2 Restriction to decorations of OCRSFs contributing to ^'dimerl^' ™) 

In this section, we characterize the restriction to decorations of OCRSFs of Qi com- 
patible with L, which contribute to ^dimer(^'^)- '^^ ^^^^ purpose we first introduce the 
following definition, see also Figure 9. 

Definition 4.1. A subset of oriented edges of the decoration x, is called an essential 
configuration compatible with L, of type cw (resp. cclwj, if it consists of: 

1. all inner edges oriented clockwise (resp. counterclockwise). 

2. one of the three following 2-edge configurations at the triangle of every non-root 
vertex Vi G i?x(L)'^- 

~{{wi,Vi),{vi,Zi)}, {{wi,Zi),{vi,Zi)}, {{vi,Wi),{wi,Zi)} in the cw case, 
~{{zi,Vi),{vi,Wi)], {{vi,Zi),{zi,Wi)], {{zi,Wi),{vi,Wi)] in the cclw case, 

3. one of the two following 1-edge configurations at the triangle of every root vertex 
Vi G i?x(L): 

— {(fWi, Vj)}, {(iWi, Zj)} in the cw case, 
—{{zi,Vi),{{zi,Wi)^ in the cclw case, 

with the additional constraint that the triangle of at least one root vertex contains 
the configuration {wi,Vi) in the cw case, and the configuration {zi,Vi) in the cclw 
case. 



z 



[w_z 




Figure 9: Essential configuration of the decoration x of type cw (left) and of type cclw 
(right). 
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When the type is not specified, we refer to the above as an essential configuration 
compatible with L, and when no confusion occurs, we omit the specification compatible 
with L. 

Proposition 12. 

1. For every decoration x, an essential configuration compatible with L is an ori- 
ented spanning forest with |i?x(L)| connected components whose roots are distinct 
elements of R^(L). 

2. Let F be an OCRSF of Qi compatible with h, which contributes to -Fdimer(^' 
then the restriction of F to every decoration x, is an essential configuration com- 
patible with L. 

Proof. Proposition 12 is a direct consequence of the next two lemmas. □ 

Lemma 13. Let F be an OCRSF of Qi compatible with L, which contributes to 
P^ijjjgj.(z, w). Then the restriction of F to every decoration x contains all inner edges 
of the decoration, one edge at triangles of root vertices, and two edges at triangles of 
non-root vertices. 

Proof. Let F be an OCRSF of Qi compatible with L. By Lemma 10, the restriction 
of F to every decoration x, is an oriented spanning forest with |i?x(L)| components, 
whose roots are distinct elements of i2x(L). This implies that: 

a. Triangles of decorations contain at most 2 edges of F. 

b. Each vertex of Qi has a unique outgoing edge of F, and each root vertex of i?x(L) 
has as outgoing edge of F the unique incident long edge (which belongs to L by 
definition of root vertices). 

c. The restriction of F to every decoration x contains 3dx — |-Rx(L)| edges. 

Let us first prove that the contribution to P^^j^^^j.{z,w) of OCRSFs having two edges at 
triangles of root vertices cancel out. Let x be a decoration of Qi, and Vi G i?x(L) be a 
root vertex of x. Then, there are three possible oriented 2-edge configurations at the 
triangle ti of the vertex Vi, which satisfy the necessary requirement of Point b. above, 
see Figure 10: 

Ti = {{wi, Vi), {zi,Vi)}, T2 = {{wi,Zi), {zi,Vi)}, Ts = {{zi,Wi), {wi,Vi)}. 

Suppose that there exists an OCRSF F^ compatible with L, having the edge config- 
uration Ti at the triangle t^. Let F^ (resp. F^) be the OCRSF F^ modified so as to 
have the edge configuration T2 (resp. T3) at the triangle ti. Then, it is straightforward 
to check that F^, F^ and F^ are OCRSFs compatible with L, with oriented non-trivial 
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T, 

Figure 10: The three possible oriented 2-edge configurations at triangles of root vertices. 



cycles in bijection (one way to prove this is to use Part 1 of Lemma 20). By Theorem 7, 
the contribution C {F'^ , F'^ , F^) of F^,F'^,F^ to -Pdimerl^'™) is: 

C{F\F\F^) := ^ ( n f^TyK^^A 11 ' ^'^^^ -^'^^^ 

Since F^,F^,F'^ have oriented non-trivial cycles in bijection, the term 0x6^^(1 ~ 
independent of j. Denote by C the contribution of edges of F^ which are 
not edges of the triangle U. Then: 



C{F\F^F^) = C n {l-z''^^^w<^^y 

{fwifvil^WijVifzifvi^ZijVi + fwifzi^^WijZifzifvi^ZijVi + fzifwiI^Zi,Wifwifvil^Wi,Vi) ■ 

Recalling that by our choice of Kasteleyn orientation, edges of the triangles are oriented 
clockwise, we have that C{F^,F^,F^) is equal to: 

c n (1 - ^'^^^«^"^^)) [fMzM-^) + (-i)/-,7^(-i)^^; + fzMwj:;;] 

= c n (1 - z^^^^w^'-^^ ) [fzJM-U +77. +7^)] 

TeF^ 

= (by definition of the vector /). 

Since this holds for every decoration x, and every root vertex Vi € i?x(L), we deduce 
that the contribution to Primer i^i"^) OCRSFs compatible with L, having two edges 
at triangles of root vertices, cancel out. 

As a consequence, if an OCRSF compatible with L contributes to -Primer ^'^^^ 
restriction of F at the decoration x contains at most, all inner edges of the decoration, 
one edge at triangles of root vertices, and two edges at triangles of non-root vertices; 
that is it contains at most dx + |-Rx(L)| -|- 2|i?x(L)'^| = 3(fx — |-Rx(L)| edges. Since 
by Point c. above, it must contain exactly 3dx — |-Rx(L)| edges, we deduce that all 
constraints must be met. □ 
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Lemma 14. A subset of edges of the decoration x, is an oriented spanning forest 
with |i2x(L)| components whose roots are distinct elements of i?x(L), satisfying the 
constraints of Lemma 13, iff it is an essential configuration compatible with L. 

Proof. By Remark 11, Lemma 14 amounts to proving that an oriented subset of the 
decoration x contains all inner edges of the decoration, one edge at triangles of root 
vertices, two edges at triangles of non-root vertices, and satisfies Condition 3', iff it is an 
essential configuration compatible with L. This is clear once we return to Definition 4.1. 

□ 

4.3.3 Essential OCRSFs of Qi 

Sections 4.3.1 and 4.3.2 above naturally yield the following definitions. 

Definition 4.2. Let r G {cw , ccwl}^ ^^^^ be an assignment of type cw or cclw to the 
vertices of G\. Define the set of essential OCRSFs of Qi compatible with L and r, to 
be the set of OCRSFs compatible with L whose restriction to every decoration x, is an 
essential configuration of type r(x). We denote this set by J^'^(^i). 

As a consequence of Lemma 10 and Proposition 12, we have the following characteri- 
zation of 

Corollary 15. A subset of oriented edges of Qi is an essential OCRSF compatible 
with L and t, iff: 

1. long edges of F are exactly those of Ij, 

2. F contains no trivial cycle consisting of short and long edges, 

3. the restriction of F to every decoration x of Qi, is an essential configuration 
compatible with L, of type r(x). 

Definition 4.3. Define the set of essential OCRSFs of Qi compatible with L, denoted 
^^(^i), by: 

and the set of essential OCRSFs ofGi, denoted J^(Gi), by: 

where C is the set of oriented edge configurations of Gi satisfying (*). 
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Corollary 16. 



i^GJ-O(Si) \e=(a;,2/)eF / TeF 



Proof. Write, J'(gi) = J^iQi) U (J^C^i))'' in the formula for P^-^^^{z,w) given by 
Theorem 7. Then use Point 2. of Proposition 12, to deduce that the contribution of 
OCRSFs of iJ^{Gi)y cancel out in P<?i^er(z, w). □ 



4.4 Statement of main result 

We can now give a precise statement of Theorem 1. 

Theorem 17. Consider a critical Ising model defined on an infinite, 7? -periodic iso- 
radial graph G. Then, one can explicitly construct essential OCRSFs of Gi counted by 
the critical dimer characteristic polynomial ^dimcr(^' from OCRSFs of Gi counted 
by the critical Laplacian characteristic polynomial PLap(z,w). 

Let us give an idea of the construction. First observe that the number of OCRSFs of 
Qi is much greater than the number of OCRSFs of Gi so that there is certainly no 
one-to-one mapping between these set of configurations. Rather, to every OCRSF of 
Gi, we assign a family of OCRSFs of Qi, which have the same reference number and the 
the same signed number of cycles, and such that the sum of the weights of OCRSFs in 
this family is equal to the weight of the original OCRSF of Gi . The family of OCRSF 
of Qi is constructed by successively adding long edges, and this construction is done 
by induction on the number of long edges added, thus allowing us to keep control over 
properties of OCRSFs in this family. To be more precise, we actually need to work on 
the graph Gi and its dual graph G^ at the same time. The operations used to add long 
edges are called licit primal/dual moves, and are valid in a context more general than 
isoradial and Fisher graphs, so that we dedicate it Section 5. The proof of Theorem 17 
is the subject of Section 6. In Section 6.1, we construct a family of CRSFs of Qi from 
each CRSF of Gi. In Theorem 29 we prove that we exactly obtain all CRSFs of Qi 
counted by the dimer characteristic polynomial; and in Section 6.2, we show that this 
construction is weight preserving. 



5 Primal/dual edge moves on pairs of dual OCRSFs 

Definitions and results of this section are valid in a context more general than isoradial 
and Fisher graphs. But, in order not to introduce too many notations, we let Gi be 
any graph embedded in the torus. 
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In Section 5.1, we first define a general OCRSF characteristic polynomial, allowing 
for weights which depend on oriented edges, and we prove a useful rewriting of this 
polynomial, using pairs of dual OCRSFs. In Section 5.3, we define and prove properties 
of licit primal/ dual edge moves, which are edge moves performed on pairs of dual 
OCRSFs, and are one of the key ingredients of the construction of Theorem 17. Licit 
primal/dual edge moves rely on a natural edge operation performed on one OCRSF 
only, which is the subject of Section 5.2. 

5.1 OCRSF chciracteristic polynomial 

We use definitions and notations introduced in Section 3.2.1. Suppose that a complex 
weight function p is assigned to oriented edges of Gi, that is every oriented edge (x, y) 
has a weight P[x,y) ^ The OCRSF characteristic polynomial, corresponding to the 
weight function p, denoted -Pocrsf(z, w), is defined by: 

Pocrsf(z,w)= p{F)]\{l-z''^'^)^N<'^)), (8) 

where p{F) = Il{x,y)eF P{x,y)- 

Remark 18. When p is a positive weight function on unoriented edges of Gi, then the 
OCRSF characteristic polynomial is simply the Laplacian characteristic polynomial. 

We now prove a useful rewriting of -Pocrsf(z, w). In order to do this, we need a few 
more facts and definitions. Denote by G* the dual graph of Gi and let us, for a moment, 
consider edge configurations as unoriented. It is a general fact that if F is a CRSF of 
Gi, then the complementary of the dual edge configuration, consisting exactly of the 
dual edges of the edges absent in F, is a CRSF of G^, with non trivial cycles parallel 
to those of F, such that primal and dual non trivial cycles alternate along the torus. 
It is referred to as the dual CRSF of F. 

Let F and F* be OCRSFs of Gi and G| respectively, then F and F* are said to be 
dual of eachother, if their unoriented versions are. This means that to a pair of dual 
CRSFs corresponds 4l'^°'^-*"^^^' ^^^^^^l pairs of dual OCRSFs. We denote by T{Gi,Gl) 
the set of pairs of dual OCRSFs, that is: 

J^{Gi,Gl) = {{F,F*)e F{Gi) X F{Gl) : F, F* are dual OCRSFs}. 

Let {F,F*) be a pair of dual OCRSFs of Gi and G^. Recall that the homology class 
H{F) of F can be rewritten as: 

H{F) = N{F)Ho{F), 

where Ho{F) = {ho{F),vo{F)) is the reference number of F, and N{F) is the signed 
number of non-trivial cycles of F defined in Section 3.2.1; and similarly for the homology 
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class of F*. Observe that F and F* have the same reference number, so that it makes 
sense to refer to Ho{F) as the reference number of the pair {F,F*). Let us denote by 
N{F, F*) the sum N{F) + N{F*). Then, we have the following rewriting of Pqcrsf- 

Lemma 19. The OCRSF characteristic polynomial of Gi can be rewritten as: 

P0CRSf(z,w) = J2 P(^^ (-z'^o(F)^.o(F))iiV(F,F*). 

(F,F*)GJ-{Gi,G*) 

Proof. Denote by Ti, • • • , T„ the OCRT components of a generic OCRSF F of Gi. By 
expanding the product of Equation (8), the polynomial -Pqcrsf (z,w) can be rewritten 
as: 

n 

Pocrsf(z,w)= Y: P(F)Il E (-z'^^^'^w'^(^'))" 

FeJ-(Gi) (£i,-,£„)e{o,i}" j=i 

Fix an OCRSF P of Gi, and let (ei, • • • , € {0, 1}". Then to F assign the following 
dual OCRSF F*: if Si = 1, take T* to be co-oriented with Ti, and if £j = 0, take T* to 
be contra-oriented. Thus, there is a one-to-one correspondence between dual OCRSFs 
of F, and sequences (si, • • • , € {0, 1}**, and: 

Si = ±^{N{Ti) + N{T*)), 

ei{h{Ti),v{Ti)) = \{h{Ti) + h{T*),v{Ti) + v{T*)). 

This implies that: 

= {^_l)\{N{Ti)+N{T*)) ^ 

ei{h{Ti),v{Ti)) = ^{N{Ti)+N{T*))iho{Ti),vo{Ti)). 

The proof is concluded by recalling that T,i=i{N{Ti), N{T*)) = {N{F),N{F*)) = 
N{F, F*) and that for every ie {!,■■■ , n}, Ho{F) = Ho{Ti). □ 

5.2 Replacing one edge of an OCRSF 

Let F be an OCRSF of Gi, and consider an oriented edge e = {x,y) of E{Gi) \ F. 
Then, since F is an OCRSF, the vertex x has a unique outgoing edge Cx of F. Define 
F{e ex} to be the edge configuration F where the edge Cx is replaced by the edge e: 

P{e,e.} :=PU{e}\{e,}. 
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In this section we describe features of the edge configuration F|g g^j., which by con- 
struction contains one outgoing edge at every vertex. 

Denote by Ti, • • • ,T„ the CRT components of F, and let 7^ be the non-trivial cycle of 
the CRT component Tj. Without loss of generality let us suppose n > 2, the case n = 1 
simply being a boundary case of the case n = 2. Then, the edge e belongs to a unique 
cylinder obtained by cutting along two neighboring non trivial cycles of F, say 71, 72, 
see Figures 11, 12, 13. Since F is an OCRSF, the vertex x (resp. y) is connected by 
an oriented edge-path p^- (rcsp. py) of F to one of the non-trivial cycles 71 or 72. Let 
F* be a dual OCRSF of F (at this point the orientation of its non-trivial cycles does 
not matter), then there is a unique non-trivial cycle 7* of F* contained in the cylinder 
^. Cutting along 7* separates the cylinder ^ into two disjoint cylinders ^1, "^2- In 
order to analyze the features of the edge configuration -P{e,ex}5 ^^^^ describe the 
edge configuration F U {e}, split according to the following cases. 



Case 1: The edge e belongs to the cylinder ^1, see Figure 11 (the case where e belongs to 
the cylinder "^2 is symmetric). Then, the vertices x and y are connected to the 
path 71, and the paths Px and py are contained in the cylinder ^1. In Case 1, 

we suppose moreover that the paths px, Py merge at a vertex m before hitting 71 
or exactly when hitting 71. Then, the edge configuration consisting of the edge e 
and the part of px and py from x (resp. y) to m consists of: 

Case la: either a trivial cycle. 

Case lb: or a non-trivial cycle parallel to 71. 




Figure 11: The edge configuration F U {e} in Case la (left) and Case lb (right). 

Case 2: As in Case 1, the edge e belongs to the cylinder '^1, but this time the paths Px,Py 
do not merge before hitting 71, see Figure 12. Then cutting along {e}U{p3;}U{py} 

separates the cylinder ^1 into two connected components, one homcomorphic to 
a disc and the other to a cylinder. We let m be the unique boundary vertex of 
the component homcomorphic to a disc, with two incoming edges. 
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Figure 12: The edge configuration F U {e} in Case 2. 

Case 3: The edge e crosses the non-trivial cycle 7*, see Figure 13. Then, up to relabeling, 
the vertex x (resp. y) is connected to the path 71 (resp. 72), and the path 
(resp. Py) is contained in the cylinder "^i (resp. '^2)- We let m be the vertex at 
which the path 7-r hits the non-trivial cycle 71 . 




Figure 13: The edge configuration F U {e} in Case 3. 

As a consequence, we obtain the following characterization of the edge configuration 

Lemma 20. Let F he an OCRSF ofG\, e = {x,y) be an edge of E{Gi) \ F, and ex he 
the unique edge of F exiting the vertex x. Then using the ahove case splitting, we have: 

1. In all cases, as long as x ^ m, then the edge configuration -^{e,ea:} o,n OCRSF 
whose oriented non-trivial cycles are in hijection with those of F, so that F and 
^{e,ex} have the same reference numher and NiF{e,ex}) — -^(-^)- 

2. When x = m, i.e. the path Px is reduced to the point m, refer to Figure 14-' 

• In Cases la and 2, a trivial cycle is created so that -F{e,e2,} not an OCRSF. 

• In Case lb: when x = m ^ 71, then a non-trivial cycle parallel to 71 is 
created, so that -P{e,ex} OCRSF with the same reference numher as F, 
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satisfying N(F^g g^j) = N{F) ± 1. When x = m G 71, then a non-trivial 
cycle parallel to 71 is created and 71 is broken up. Thus, -F{e,ea;} is an OCRSF 
with the same reference number as F, satisfying N{F^gg^y) = N{F) ±2 or 

iV(F{,,,^}) = iV(F). 

• In Case 3; as long as the number n of non-trivial cycles is > 2, then the 
non-trivial cycle 71 is broken up, so that -^{e,ea:} '^f^ OCRSF with the same 
reference number as F, satisfying N{F^g ,,^y) = N{F)zizl. When n = 1, then 
the unique non-trivial cycle is broken up, and a non-trivial cycle orthogonal 
to 71 is created. Thus -^{e,ex} ^ different reference number than F. 




Case 1 a Case 1 b Case 2 Case 3 



Figure 14: The edge configuration F|g g^j, when x = m. 
5.3 Licit primal/dual edge moves 

Let (F, F*) be a pair of dual OCRSFs of Gi and G\. We consider an edge ei = (.xi, yi) 
of E{Gi) \ F, and an edge 63 = {x2,y2) of E{Gl) \ F*. Then, ei is the dual of an 
oriented edge el = {xl,yl) of F* , and 62 is the dual of an oriented edge 62 = {x2,y2) 
of F. 

Definition 5.1. We say that the pair of primal and dual oriented edge configurations 
(-^{61,62}' -^{e* e*}) obtained from (F, F*) by a primal/dual edge move. This move is 
called licit, whenever: 

xi = X2 and xl = X2. (9) 

Proposition 21. The pair (Fjg^ ejji ^{e* e*}) consists of dual OCRSFs of Gi and G^, 
iff the primal/ dual edge move is licit. When this is the case, 

1. either the pair (F[g^ g^j., F^^. g.^) has the same reference number as {F,F*), and 

N{F{e,,e,},F{e*,e*,}) = N{F,F*), 

2. or the pairs (Ffgj^g^}., F|'g, ^,^) and {F,F*) have different reference numbers, and 

NiF^e,M^Fl^.^^.y) = N{F,F*) = 0. 
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Proof. Let us first show that Condition (9) is necessary. By construction of F^ei,e2}' 
the outgoing edge ei is added at the vertex xi. Since F is an OCRSF, the vertex xi 
has a unique outgoing edge of F. Suppose that xi X2, this impHes that 62 7^ Sx-^ 
so that the vertex xi has two outgoing edges in the edge configuration -^{ei,e2}- 
Remark 2, this forbids -F{ei,e2} from being an OCRSF. A similar argument holds for 

Let us now suppose that xi = X2 and xl = X2, see Figure 15. 



,,* e- 



Figure 15: The edges ei, 62, e^, e^, when xi = X2, = a;|. 

In the sequel, we use the notations of Section 5.2. Suppose that xi 7^ mi, and .T2 7^ 
Then, by Point 1. of Lemma 20, the edge configuration ^^{61,62} (resp. F^^t, ) has the 
same reference number as F (resp. F*) and -^(-^{ei,e2}) ~ ^(F) (resp. N{F^^, = 
N{F*)), so that Point 1. of Proposition 21 is clearly satisfied. 

Suppose now that xi = mi, then in Cases la and 2, the edge-path py-^ hits the vertex 
xi, and {ei} Li {py^} contains a trivial cycle, see Figure 16. In the dual graph, the edge- 
path Px* must enter this trivial cycle, which is impossible since Px* must also connect 
xl to one of the dual non-trivial cycles, so that this case cannot occur. 




i 

y2 



Figure 16: In Cases la and 2, when xi = mi, the edge-path hits xi, and {ei}U{py^} 
contains a trivial cycle. 

In Case lb, the edge-path py^ hits the vertex xi, and {ei} U {py^} contains a non- 
trivial cycle parallel to 71, sec Figure 17 (left). Suppose now that xi = mi ^ 71, 
then in the dual graph, the edge path py* must hit X2 (implying that X2 = 1^2), and 
{62} U {Py*} must contain a non-trivial cycle parallel to 71, implying that the non- 
trivial cycles contained in {ei} U {py^} and {62} U {py*} are parallel and in opposite 
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directions, sec Figure 17 (right). Thus, by Point 2., Case lb of Lemma 20, we know that 
(F{ej e2}) -^{e* e*}) ^^^^ OCRSFs of Gi and Gl, with the same reference number as 
{F,F*), such that: 




Figure 17: Left: In Cases 16, when xi = mi, the edge-path pj^j hits xi, and {ei}U{py^} 
contains a non-trivial cycle. Right: when xi = mi ^ 71, the edge path py* must hit X2, 
and {62} U {Py*} contains a non-trivial cycle, with opposite direction. 

In Case 16, when xi = nii G 71, then {ei} U {py,^} and 71 must be co-oriented, so that 
by Point 2., Case 16 of Lemma 20, we know that Fjg^ g^} is an OCRSF, in which the 
non-trivial cycle 71 is broken up and the non-trivial cycle {ei}U {py^ } is added, implying 
that -P{ei,e2} has the same reference number as F, and iV(F{ej g^}) = N{F). In the dual 
graph, wc have Xo m%, so that by Point 1. of Lemma 20, wc know that , „,i is 
an OCRSF whose non-trivial cycles arc in bijcction with those of F. Summing the two 
contributions -^(-^{ei,e2}) N{F^^, we deduce that Point 1. of Proposition 21 is 
satisfied. 

In Case 3, when xi = mi, the edge must belong to 7* (implying that xl = X2 = mg), 
and the paths 7* and 71 must be in opposite directions, see Figure 18. If the number 
of non-trivial cycles of F is > 2, then by Point 2., Case 3 of Lemma 20, we know that 
-^{61,62} (I'esp. -F|g, is an OCRSF, in which the non-trivial cycle 71 (rcsp. 7*) is 
broken up. As a consequence, (-^{61,62}) -^{e* e*}) dual OCRSFs of Gi and Gl with 
the same reference number as [F, F*) such that: 

^(^{ei,e2},^{e*,e^}) = ^(^) ± ^ + ^(^*) ^ 1 = N{F,F*). 

When the number of non-trivial cycles of F is 1, since 7* and 71 must be in op- 
posite directions, we know that N{F,F*) = 0. Moreover, by Point 2., Case 3 of 
Lemma 20, we know that -F{ei,e2} (resp. F^^, ^ty) is an OCRSF, in which the non- 
trivial cycle 71 (resp. 7*) is broken up and in which a non-trivial cycle orthogonal 
to the original one is created. Since these must have opposite directions, we deduce 
that the pair (-P{ei,e2}' -^{V e*}) ^ different reference number than {F,F*), and 
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Figure 18: In Case 3, when xi = mi, the edge el must belong to 7*, and the paths 7* 
and 71 must be in opposite directions 

-^(-^{ei,e2}>-^{e* e*}) = = N{F,F*), implying that Point 2. of Proposition 21 is satis- 
fied. ' " ' □ 

6 Proof of Theorem 17 

Let G be an infinite, Z^-periodic isoradial graph, Q be the corresponding Fisher graph, 
and Gi, Qi be their respective fundamental domains. As a consequence of Lemma 19 
and using notations introduced in Section 4.1, we deduce that the critical Laplacian 
characteristic polynomial PLap(z,w), and the critical dimer characteristic polynomial 
P^jjjjgj.(z, w), can be rewritten as: 



In Section 6.1, we prove the first part of Theorem 17 by constructing pairs of essential 
OCRSFs of Gi and counted by P^^^^^{z,^N) from pairs of dual OCRSFs of d and 

Gl counted by PLap(z,w), sec Theorem 29. Then, in Section 6.2, wc prove that this 
construction is weight preserving, thus ending the proof of Theorem 17. 

6.1 Explicit construction 

Let us start by giving the general idea of the construction. Let (F,F*) be a pair of 
dual OCRSFs of Gi and GJ, then to (F,F*) we assign a family (S(F, F*) consisting of 
pairs of dual essential OCRSFs of Qi and Ql, such that: 




1- U(F,F*)e^(G„G^) «5(F, F*) = J^iGuGt). 
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2. When (Fi, FJ) ^ (F2, F^), then 5(Fi, FJ) n S{F2, F^) = 0. 

3. For every {F,F*) G 5(F,F*), then 

• either the reference number of {F,F*) in Qi, Q\ is equal to the reference 
number of (F,F*) in d, G\, and Ng^^g*{F,F*) = Ng,,gi{F,F*), 

• or the pairs {F,F*) and (F,F*) have different reference numbers, and 

Ng,,g^^{F,F*) = Ng,,gi{F,F*) = 0, 

where we have added a subscript to the signed number of cycles N(-) to indicate 
on which graph it is computed. 

For every pair (F,F*) of dual OCRSFs of Gi and G*, loosely stated, the family 5(F, F*) 
is constructed as follows. Let ei, • • • ,em, be an arbitrary labeling of unoriented edges 
of E{Gi) \ F. For A; G {0, • • • ,m}, we let Jfc = {(ii, • • • ,ik) G {!,••• ,m}'' | 1 < ii < 
■ ■ ■ < h ^ with the convention that Jfc = 0, when A; = 0. Then, 

m 

5(F,F*)=|J IJ JF(F,F*),{eH,-,eiJ(g^^g*)^ 
A;=0 (n,---,ifc)eJfc 

where F^^'^ '^^k^(^Q^^ Q*'j is constructed by induction on k using licit primal/dual 

moves, introduced in Section 5. In Section 6.1.1, we prove all results needed for the 
initial step of the induction. Then, in Section 6.1.2, we specify licit primal/dual moves 
to pairs of essential OCRSFs of Qi and Q*. This allows us to precisely define the set 
jr(F,F*),{e,^,--,e,J^g^^ g*^ in Section 6.1.3. Finally, we state and prove Theorem 29 
establishing that we exactly obtain all pairs of dual essential CRSFs of Qi and Q^. 

6.1.1 Initial step of the induction 

Let F be a subset of oriented edges of Gi such that each vertex has exactly one outgoing 
edge of F. Considering F as a subset of long edges of Qi defines, for every decoration x, 
a unique root vertex vi(x). From now on, whenever no confusion occurs, we omit the 
argument x. For the whole of this section, we let r G {cw, cdw}^'^'^^^ be an assignment 
of type cw or cclw to vertices of Gi . 

Lemma 22. A subset of oriented edges F ofQi is an essential OCRSF compatible with 
F and t, iff: 

1. F is an OCRSF of Gi, 

2. long edges of F are exactly those ofF, 

3. the restriction of F to every decoration x contains: 

• all inner edges with the orientation induced by the type t(x), 
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• the edge {w\,vi) at the triangle ti if the decoration is of type cw, or the edge 
{z\,v\) if it is of type cclw, 

• for every i 1, any of the three possible 2-edge configuration of Defini- 
tion 4-1, at the triangle ti, with the orientation induced by the type, see also 
Figure 9. 

Moreover, when this is the case, oriented non-trivial cycles of F in Qi are in bijection 
with oriented non-trivial cycles of F in Gi. 

Proof. Since each vertex of Gi has a unique outgoing edge of F, we know that F is an 
OCRSF of Gi, iff it contains no trivial cycle. Moreover, we know that F is an essential 
OCRSF of Qi compatible with F and r, iff it satisfies Conditions 1.2.3. of Corollary 15. 
Since every decoration of Qi has a unique root vertex, Condition 3. can be rewritten 
as Condition 3. of Lemma 22 above. Thus, it remains to show that F contains no 
trivial cycle consisting of long and short edges of Qi, iff F contains no trivial cycle of 
Gi, or equivalently F contains a trivial cycle consisting of long and short edges iff F 
contains a trivial cycle. Because of the geometry of Qi, the direction from left to right is 
straightforward. Conversely, suppose that F contains a trivial cycle. Then, since every 
vertex of Gi has a unique outgoing edge of F this implies that the trivial cycle must 
be co-oriented. Now, consider an oriented configuration of Qi satisfying Conditions 
2. and 3. above. Since every decoration of Qi has a unique root vertex, this means 
that every vertex of the decoration x is connected by a unique path to the root vertex 
wi(x), in particular this holds for every vertex of type 'u' of the decoration. Thus, if 
F contains a trivial cycles, then F contains a trivial cycle consisting of long and short 
edges. A similar argument shows that if F is an OCRSF of Qi, its non-trivial cycles 
are in bijection with those of F in Gi. □ 

Figure 19 below illustrates in an example all essential OCRSFs of Qi compatible with 
F and r. 




Figure 19: Left: an OCRSF F of Gi and an assignment of types to vertices of Gi. 
Right: essential OCRSFs of Qi compatible with F and r. 
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We now study properties of dual OCRSFs. Observe that the dual graph GJ is a 
subgraph of Q*, see also Figure 20 below. Since inner edges of the decorations are 
always present in essential OCRSFs of Qi, we omit their dual edges in our picture 
of 
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Figure 20: Left: the graph Gi (dotted lines) and it dual G\ (full lines). Right: the 
graph Qi (dotted lines) and its dual Ql (full lines), without dual edges of inner edges 
of the decorations. 

Let F be an OCRSF of Gi and F be an essential OCRSF of Qi compatible with F 
and T. Suppose for the moment that edges are unoriented, and denote by F* the dual 
CRSF of F m.gi, and by F* the dual CRSF of F in G^ Then, the restriction of F* 
to dual long edges is F*, and since G\ is a subgraph of Q^, the CRSF F* is a subgraph 
of F*. 

Moreover, by Lemma 22, we know that non-trivial cycles of F in Qi and F in Gi are 
in bijection, implying that non-trivial cycles of F* in QI and F* in G\ are also in 
bijection. As a consequence, the non-trivial cycles of F* in Ql are exactly the non 
trivial cycles of F* in G\, and branches of F* in G\ are also branches of F* in Q^- 
Recalling that OCRSFs are obtained from a CRSF by orienting branches towards the 
non-trivial cycles, and orienting each of the non-trivial cycles in one of the two possible 
ways, we have shown the following lemma describing oriented versions of F* and F*. 

Lemma 23. Let F he an OCRSF of Gi and F he an essential OCRSF of Qi compatihle 
with F and r. Then, if F* is a dual OCRSF of F, the restriction F* of F* to dual 
long edges is an OCRSF of G\, and oriented non-trivial cycles of F* in Ql are exactly 
those ofF* in G^. Conversely, let F* he a dual OCRSF ofF, then there is a unique 
dual OCRSF F* of F in Ql such that the restriction of F* to dual long edges is F*. 

Figure 21 below illustrates a dual OCRSF F* of Figure 19, and the unique dual OCRSF 
of F of Figure 19, whose restriction to dual long edges is F*. 



34 



Figure 21: Left: A dual OCRSF F* of F of Figure 19. Right: the unique dual OCRSF 
F* of F of Figure 19 whose restriction to dual long edges is F*. 

Definition 6.1. Let (F,F*) be a pair of dual OCRSFs of d and Gl- Then, define 
pr,(F,F*)^g^ g*) to be the set of pairs {F,F*) of essential dual OCRSFs ofQi and Gf, 

com,pa,tihle with (F, F*) and r, that is the set of pairs (F, F*), such that F is an essential 
OCRSF of Qi compatible with F and r, and F* is the unique dual OCRSF of F whose 
restriction to dual long edges is F*. 

As a consequence of Lemmas 22 and 23, we have the following. 

Corollary 24. The set F^^^^''^*\gi,gi) is well defined and non-empty. Moreover, the 
reference number in g\ of every pair {F, F* ) in this set is equal to the reference number 
in Gi of (F,F*) in Gi, and 

Ng,,g*iF,F*)=NG,,GliF,F*). (10) 

6.1.2 Essential moves and reverse moves 

Let L be a subset of oriented edges of Gi also considered as a subset of long edges 
of gi, satisfying (*) of Section 4.3.1 (that is every decoration of has at least one 
root vertex of L). For the whole of this section, we let r G {cw,cdw}^^^^^ be an 
assignment of types to vertices of Gi. Suppose that there exists an essential OCRSF 
F of gi compatible with L and r, and let F* be a dual OCRSF of F. 

In this section, we characterize licit primal/dual moves performed on {F,F*), yielding 
a pair of essential dual OCRSFs of and gf, such that the first component either has 
an additional long edge {essential move) or a long edge of L removed {essential reverse 
move). 

Notation. For the remainder of the paper, we need to introduce a specific notation 
for oriented edges of Gi . Fix an arbitrary orientation of edges of Gi , and denote by e a 
generic unoriented edge, by +e the oriented edge compatible with the fixed orientation, 
and by — e the reverse oriented edge. 
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Essential moves 



Assume that L 7^ E(Gi), and let e be an unoriented edge of E{Gi) \ L. We now 
characterize licit primal/dual moves performed on (F, F*), which yield a pair of essential 
dual OCRSFs of Qi and such that the first component is an essential OCRSF 
compatible with L U {±e} and r, i.e. belongs to F'^'^^^^^^ {Qi). 

Suppose first that we want to add the edge +e to F, corresponding to a long edge 
{vk{^)-,Vi{y)) of Qi, for some G {1, • • • , dx}, ^ £ {1, • ' ' > dy}, see Figure 22. Assume 
for the moment that x has type r(x) = cw, and let us omit the arguments x and y, 
recalling that the index k refers to x, and the index I to y. 

Since {vk,V£) is absent in F, its dual edge {x*,y*) is present in F*. Let us first handle 
Case A, where x* is to the right of {vk,vi), see Figure 22 (first two columns). Since 
F is an essential OCRSF it has one of the three possible 2-edge configurations at the 
triangle of the vertex Vk, with the orientation induced by the type r(x) = cw. Since 
we want the resulting configuration to be an essential OCRSF compatible with r, this 
implies that after the move, the triangle must contain either the edge {w^, v^) or the 
edge {wk,Zk)- Then, see Figure 22 (first column), our only choice is to remove the edge 
{vk-,Zk) in Cases AI and All, and the edge {vk,Wk) in Case AIII. This constraint fully 
determines the dual edges which are removed/added, and we only allow the move when 
it is licit, i.e. when edges involved satisfy Condition (9), that is in Cases AI and All, 
see Figure 22 (second column). In Case B, when x* is to the left of {vk,Vi), a similar 
argument holds and yields the last two columns of Figure 22. 




Figure 22: Adding the long edge {vk,Vi). First two columns: x* is to the right of 
{vk,Vi). Last two columns: x* is to the left of {vk,ve). 

By symmetry, when x has type r(x) = cclw, if x* is to the right (resp. to the left) of 
{vk,V£), we get Case B (resp. Case A) with 'w's and 'z's exchanged. Using symmetries 
again, the case where we add the edge {vi,Vk) to the OCRSF F is also handled by 
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Cases A and B. It is important to note that in each case at most one primal/dual move 
is allowed. 

Definition 6.2. When the above move is allowed, it is called an essential move, and 
I = .1., F| denotes the resulting pair of OCRSFs. When it is not 

allowed, we set by convention, {F,F*)^±.^ .y = 0. 

Then, Proposition 21 and the construction of essential moves immediately yields the 
following. 

Proposition 25. Let F be an essential OCRSF of Qi compatible with L and t, and 
F* be a dual OCRSF of F. Then, the pair (F, F*){±e, } consists of: 

• the empty set when the essential move is not allowed, 

• when the essential move is allowed, a pair of dual essential OCRSFs of Qi and 
Ql, such that F is an essential OCRSF compatible with L U {±e} and t, and 

1. either the pair -^*){±e, } same reference number as {F,F*), and 

N{{F,F*){^,^.y) = N{F,F*), 

2. or the pairs -^^*){±e,-} (^i^d {F,F*) have different reference numbers, and 

N{{F,F*){±,^.y) = N{F,F*) = 0. 

Moreover, the orientation of edges not involved in the move remains unchanged. 
Essential reverse moves 

Assume now that L 7^ 0, and let ee = (x, y) be an oriented edge of L, where e G {— , +}. 

We now characterize licit primal/dual moves performed on {F,F*), which yield a pair 
of essential OCRSFs of Qi and QI, such that the first component is an essential OCRSF 
compatible with L \ {ee} and r, i.e. belongs to J"^'^^^^'^^ (Qi) . 

The edge ee corresponds to a long edge {vk{x.) , ve{y)) of Qi, see Figure 23. Assume 
for the moment that x is of type r(x) = cw, and let us omit the arguments x and y 
onwards. Since F is an essential OCRSF, it contains either the edge {wk,vi^) or the 
edge {wk,Zk) at the triangle tk- Note that the edge ZkVk is always absent, so that its 
dual edge is present in F*. Let us first handle Case C, where this dual edge is oriented 
towards the triangle t^. 

Since we want the resulting configuration to be an essential OCRSF compatible with r, 

this implies that after the move, the triangle t^ must contain one of the three possible 
2-edge configurations, with the orientation induced by t(x) = cw. Then, see Figure 
23, when the triangle t^ contains the edge {wk,Vk), our only choice is to add the edge 
{vk,Zk) (Case CI). When it contains the edge {wk,Zk), we can either add the edge 
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Case C Case D 



Figure 23: Removing the long edge (vk,Vi) when the dual of the edge ZkVk is oriented 
towards (away from) the triangle tk in Case C (Case D). 

ivk-,Zk) (Case CII) or the edge {vk,Wk) (Case CIII). We only allow this move when it 
is licit, that is in Cases CI and CII. 

In Case D, where the dual of the edge z^Vk is oriented away from the triangle t^, a 
similar argument holds, and yields the last two columns of Figure 23, the move is only 
allowed in Case DHL It is important to note that in each case, at most one move is 
allowed. By symmetry, if x has type cclw, Cases C and D hold, with ^w^s and 'z's 
exchanged. 

Definition 6.3. When the above move is allowed, we refer to it as an essential reverse 
move, and denote by (F, = (F^.^^e}; -^{ee*,-}) the resulting pair of OCRSFs. 

When the move is not allowed, we set by convention (F, F*)|. £e} = 0- 

The analogous of Proposition 25 also holds for essential reverse moves, and we do not 
write it out explicitly. The next lemma relates essential moves and reverse ones. 

Lemma 26. Let (F, F*) be a pair of dual essential OCRSFs of Gi and Q\, such that 
F is compatible with L and r. 

1. Let e be an unoriented edge of E{G\) \ L, and e G {— , +}. Then, if an essential 

m,ove can he performed on {F,F)*, yielding a pair {F, F*)^^^ .^, the essential re- 
verse primal/dual move, which removes the long edge ee from {F, F*)^ge .y, can 
also be performed and: 

((F,F*){,e,.}){,ee} = (i^,^*)- 

2. Let ee be an oriented edge ofh, where e G {—,+}• Then, if an essential reverse 
primal/dual move can be performed on {F,F*), yielding a pair {F, F*)^. ,,^y, the 
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essential primal/dual move, which adds the long edge ee to (F, ^gj., can also 
performed, and: 

((F,F*){.,,e}){ee,.} = (^'^*)- 

Proof. From Figures 22 and 23, we immediately check that the moves AI/CI, AII/CII, 
BIII/DIII are inverse of eachother. □ 

6.1.3 Explicit construction 

Let (F,F*) be a pair of dual OCRSFs of Gi and G^, and let ei.-- - ,e„i, be an 
arbitrary labeling of unoriented edges of E{Gi) \ F. We let Jk — {(ii,--- G 
{!,••• , Txi\^ 1 1 < n < • • ■ < 'ifc < m)-, with the convention that = 0, when A; = 0. 
Let r G {cw, cclw)^ ^'^'^^ be an assignment of types to vertices of G\. 

Wc now define by induction, for every A; G {0, • • • , m} and every (ii, • • • , ife) G Jfc, the 
set, 

jT-,(F,F*),{eij ,- ,eij^ ;•} (g^ ^ g*) . 

Initial step : /"^'(^'^'^'^gi, ^i*) = {Qx.QD , defined in Section 6.LL 

Induction step: for _7 = 1, • • • , fc, 

U {(F,F*){+e,^.,.}U(F,F*){_e^^.,}}, 
where (-^^, -F*){±ei.,-} is defined in Section 6.1.2. Then, we have: 

Proposition 27. The sei J^'^^'^*^''^^'i''"'^'fc''^(^i,^*) is non empty, independent of the 
order in which edges are added, and consists of distinct pairs of essential dual OCRSFs 
of Gi and Ql, compatible with r. Moreover, for every pair {F,F*) in this set, 

1. either the reference number of {F,F*) in Qi is equal to the reference number of 
(F,F)* in Gi, and: 

Ng,,g*^{F,F*) = NG„Gt{F,F*), 

2. or the reference numbers differ, and 

Ng„g^^{F,F*) = Ng,,gi{F,¥*) 
Remark 28. Note that the set ■'^*fc''i^(^i,^j|") can be rewritten as: 
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where J^'^'^^'^*^'i^n^n''" '^^k^^k''} (^g^^g*') is defined by induction as follows. The initial 
set is ^if ), and for every j e {!,■■■ , k}, 

U {(i^,i^*){.,e,.,.}}- 

(F,F*)eJF"'''''''*^'^'^i''*i''''*i-i%-i'\gi,ej) 

Proof. As a consequence of Remark 28, it suffices to show that for every (e^j , • • • , 6,^,) G 
{-,+}'=, the set 7^'(^'^'')'K^n'-'%''ifc'-}(^^,^*) satisfies Proposition 27. Let us first 
prove that the initial set J-"^'^^'^*\Qi, Q^) satisfies all statements. By Lemma 22 and 23, 
it is clearly non empty, and consists of pairs of essential dual OCRSFs of and g^, 
compatible with r; and by Corollary 24, it satisfies Point 1. 

Returning to the definition of essential moves and to the characterization of OCRSFs of 
^r,(F,F*)(^^^^*)^ we deduce that the set J^'C^'^*^'^'^^^'- '^^k^ik'-}{g^,gi) is also non- 
empty. By Proposition 25, we know that it consists of pairs of essential dual OCRSFs 
of gi and gi compatible with r. Moreover, since the orientation of edges not involved 
in the move remains unchanged, we deduce that the induction is in fact independent 
of the order in which edges are added. Proposition 25 also implies that either Point 1. 
or 2. is satisfied. □ 

Definition 6.4. 

fc=0 {{e,^ ,ei^ }:(«!,■■■ ,ik)&Jk} 

s^^^^'\guGl)= U s-'(^^^'\g,,gi). 

The next theorem proves that we exactly obtain all pairs of dual CRSFs of and gf. 
Theorem 29. 

2. When (Fi,F|) + (F2,F^), then ^^^I'^^^l^i, ^i*) n ^(f^.f*)^^^^ ^ 0_ 

Proof. The inclusion C of Point 1. is a direct consequence of Proposition 27. 

Consider a pair {F,F*) of dual essential OCRSFs of and g*, and let us show 
that there exists a pair of dual OCRSFs (F,F*) of Gi and G\, such that {F,F*) G 
S^^'^*\gi, gi). By definition of essential OCRSFs of ^i, there exists a subset of edges 
L of C, and r G {cw, cclw}^^'^^\ such that F is an essential OCRSF compatible with L 
and r, and F* is a dual OCRSF of F. Recall that the set L defines, for every decoration 
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X, a set of root vertices i?x(L). Let us fix a decoration x, assume tliat r(x) = cw {the 
case where r(x) = ccZ?« being similar), and omit the argument x in the sequel. By 
Corollary 15, the restriction of F to the decoration x consists of: 

- all inner edges oriented clockwise, 

- one of the three possible 2-edge configuration at the triangle of every non-root 
vertex, with the appropriate orientation, 

- one of the two following 1-edge configurations at the triangle of every root vertex 
Vi e -Rx(L): 

{{Wi,Vi)},{{Wi,Zi)}, 

with the additional constraint that the triangle of at least one root vertex contains 
the configuration (wi,Vi). 

Let us now study properties of the dual OCRSF F* at the decoration x, sec also 
Figure 24. Denote by c* the dual vertex at the center of the decoration, and by 
tl, - ■ ■ , the dual vertices at the center of the triangles ti, ■ ■ ■ ,t(i^. As inner edges 
are always present in the OCRSF F, we omit their dual edges in our representation of 
the dual graph. 




Figure 24: The restriction of F and F* to a decoration. 

Since F* is an OCRSF, there is exactly one edge {c*,t*) exiting the vertex c*, we set 
by convention i = 1. Then, since there must also be one edge exiting the vertex t^, this 
means that F must contain the edge (wi,vi), implying that vi is a root vertex of L, 
and that the long edge eiei, where ei G {—,+}, whose initial vertex is vi, belongs to L. 
Repeating this for every decoration of Qi defines a subset of oriented edges F of Gi, 

F= (J {ei(x)ei(x)}, (11) 
xeV(Gi) 

such that every vertex of Gi has exactly one outgoing edge of this subset. 

Suppose that the decoration x has another root vertex Vi (i ^ 1), and let eje^ be the 
edge of L whose initial vertex is Vj. Then, at the triangle ti, the OCRSF F consists of: 
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• either the edge {wi,Zi); then, the dual of the edges WiVi, ZiVi belong to F*. A 
priori, there are two possible orientations for the dual edges, and which one it is 
is fixed by F*. In each of the two cases, either the essential reverse move CII or 
Dili can be performed on the pair (F,F*), yielding a pair (F, F*)^.^^.^.j of dual 
essential OCRSFs of Qi and Q^. 

• or the edge {wi, Vi); then, the dual edge (t*, c*) belongs to F* and must be oriented 
towards c*. Indeed, otherwise the vertex c* would have two exiting edges of F* 
which contradicts the fact of being an OCRSF. This implies that the dual of the 
edge ZiVi, which also belongs to F* is oriented towards t*. As a consequence, the 
essential reverse move CI can be performed on the pair (F,F*), yielding a pair 
(F,F*){.^,.e,} of dual essential OCRSFs of and Gl 

Since the orientation of edges not involved in the move remains unchanged by essential 
reverse moves, we can repeat this for every root vertex different from vi of the decoration 
X, and for every decoration x of ^i. By the analogous of Proposition 25 for essential 
reverse moves, this yields a pair {F,F*) of dual essential OCRSFs of Qi and Q'l, such 
that F is compatible with F and r. Then, by Lemma 22, this implies that F is an 
OCRSF of Gi. Let F* be the restriction to dual long edges of F*. Then, by Lemma 23, 
F* is a dual OCRSF of F, and we deduce that: 

{F,F*)eJ^'^^'^''\gi,gi). 

By Lemma 26, if an essential reverse move can be performed on a pair of dual OCRSFs, 
removing an oriented long edge, then the essential move adding the same long edge can 
be performed to recover the original pair. Applying this recursively, we deduce that: 

{F,F*) G j^'(^'^'^'i^\^y{g,,gi) c 

thus proving Point 1. of Theorem 29. 

Let us now prove Point 2. Suppose that there are two distinct pairs (Fi, FJ) and (F2, F2) 
of dual OCRSFs of d and G\, such that: 

{F,F*) e s^^'^^*^\gi,gi)ns^^^'^*^\gi,gi). 

Let us return to properties of {F, F*) at a fixed decoration x of ^1. In F*, the dual edge 
of the edge zivi is oriented away from t^. Now, suppose that the edge eiei is added by 
an essential move. This implies that in the dual graph, the dual of the edge zivi must 
be added. Referring to Figure 22, which describes possible essential moves, we see that 
this dual edge is then oriented towards ti(x), which is a contradiction. Thus, eiei must 
be an edge of the original OCRSF. Since essential moves do not change the orientation 
of edges not involved in the move, we repeat this argument for every decoration, and 
deduce that: 

Fi = F2 = F, 
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where F is defined in Equation (11). As a consequence (F, F*) € J="''^^'^*^'^^\^\gi,gi)n 
jrr,{F,F*),{L\F}(g.^ g.*)^ where and F* are two distinct dual OCRSFs of F. This 
means that there exists € ^i*) and {F2,F^) € -F^'^^'^^H^i, 

such that {F, F*) is obtained from each of (Fj, Fi)* by successively adding the same set 
of oriented long edges L\F with essential moves. By Part 1. of Lemma 27, if an essential 
move can be performed to add an oriented long edge, then the reverse move can also 
be performed to recover the original pair. Moreover, looking at Figure 9, describing 
possible essential reverse moves, we see that there is at most one essential reverse move 
for removing a given oriented long edge. Thus, the same essential reverse moves are 
performed on {F, F*) to recover (Fi, Fj*) and (Fj, F2*), thus implying that this is in fact 
the same pair. Since F* is the restriction to dual long edges of F* (i = 1, 2), we deduce 
that F| = F^. □ 



6.2 The construction is weight preserving 

In this section, we state and prove Theorem 30, which establishes that the construction 
is weight preserving: we show that the sum of the weights of all pairs of essential dual 
OCRSFs of ^1 and in S^'^^^*\gi,g*) is equal, up to a constant which only depends 
on the graph d, to the weight of the pair of dual OCRSFs (F,F*) of Gi and G^. 
As a consequence we recover, by an explicit computation. Theorem 2 of [BdTlOb], see 
Corollary 31. Note that the constant could not be explicited in the proof of Theorem 
2, but could only be recovered a posteriori in [BdTlOa]. 

Theorem 30. Let {¥,¥*) be a pair of dual OCRSFs of d and GJ. Then, 
(F,F*)e5(F.F*)(e;i,6i*) \{x,y)eF J 

= C[ n tan^xJ (-z'^°(P)w"o(P))l^(^'P*), 

\(x,y)eF / 

where C = 24l^(Gi)l+lV(Gi)| nxye£;(Gi) sm2(^) cos^^y. 

Proof. Using notations introduced in Section 6.1.3, let us recall that: 

m 

s(^'^'\g,,gi) = U U U .F^'(^'^*)'{^H.-.e,j(g^,g*). 

By Proposition 27, we know that for every pair (F, F*) of dual essential OCRSFs of 

s(-^'^'\gugty. 
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so that we only need to handle the weights which do not involve the coefficients z 
and w. Consider r G {cw,cdw}^^'^^\ a fixed assignment of types to vertices of Gi, 
and let us denote by Wk the weighted sum (excluding z and w) of configurations of 

jrr,(F,F*),{eij,-,eij^}^g^^ -yy-g compute Wfc by induction on k. 

In all computations below, we use the definition of the Kastclcyn orientation of Sec- 
tion 3.1, the definition of the function {fx)xeV{gi) of the rhombus half-angles in 
M/47rZ of Section 4.2.1. 

Computation of Wp . Recall that (^i, Q*) consists of all pairs {F, F*) of dual 

OCRSFs of Qi and Q^, such that F is an essential OCRSF compatible with F and r, 
characterized in Lemma 22, and F* is the unique dual OCRSF whose restriction to 
dual long edges is F*. Recall also that, for every decoration x of ^i, fi(x) denotes the 
unique root vertex induced by F. 

1. Contribution of long edges of F. Let (x, y) be an oriented edge of F, and let 
{vk{x),Vi{y)) be the corresponding long edge of Qi. Denote by ^xy the rhombus 
half-angle of xy. Then the contribution of this edge to Wq is: 

fvk fvi ,Vl ■ 

By definition of the Kasteleyn matrix and of the vectors /, we have: 

K -£ cot ^ 

= e ' 2 - e * 2 , /^^ = e * 2 - e * 2 . 
Moreover, by definition of the rhombus half-angles in R/47rZ: 

fve = i^Vk,vtfvk-i 
SO that fy^ = —i£vk,vifvk- 3, consequence: 

fvJ7eKvk,ve = -i cot ^ 1/^^ P = -i cot ^ 2(1 - cos ^^y) 

= -2isin6'xy. (12) 

2. Contribution of inner edges of decorations. Let x be a decoration of Qi, and let 

■Wj{x)zj^i{x.) be a a generic inner edge of x. 

o If t(x) = cw, then the edge is oriented from Zj+i to Wj, and its contribution to 
Wo is: 

£zj + l,Wjfzj + lfwj- 

Moreover, by definition of the angles in R/AttZ, we have fzj+i = —£wj,zj+ifwj, so 
that the contribution is: 

~£zj + l,Wj£wj,Zj + lfwjfwj — 1- 
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o If r(x) = cclw, then the edge is oriented from wj to ^j+i, and its contribution 
to Wo is: 

3. Contribution of triangles of root vertices. Let x be a decoration of Qi, and ti(x) 
be the triangle of the root vertex fi(x). Denote by ^i(x) the rhombus half-angle 
of the long edge incident to Vi (x) . 

o If t(x) = cw, then the triangle ti contains the edge {wi,vi), so that its contri- 
bution to Wo is: 

fwJviK^uv^=e * 2 (e' 2 -e* 2 ) = l-e = 2z sm(^— j e *2. 

o If r(x) = cclw, then the triangle contains the edge {zi,vi), so that its 
contribution to Wo is: 

fzJv,K,,,v, = e-^— (e'— - e*— ) = -1 + e^^^ = 2^sin(^je^^. 

4. Contribution of triangles of non-root vertices. Let x be a decoration of Qi, tj{x) 
(j 7^ 1) be a triangle of a non-root vertex, and 6j be the rhombus half-angle of 
the long edge incident to Vj(x). Then, tj(x) contains any of the three possible 
2-edge configurations, with the orientation induced by the type. Thus, 

o If r(x) = cw, its contribution to Wo is, see also Figure 9: 



fwj fvj -^Wj ,Vj fvj fzj Kvj,Zj ~\~ fvj fzj -^Vj ,Zj fwj fzj -^Wj ,Zj ~^ fvj fwj -^Vj ,Wj fwj fzj -^Wj ,z 
~ fwjfzjfvjifvj fzj ~l" fwj) 

= fwj fzj fvj {2fwj ) (by definition of the vector /) 
= VzJv, = -2e'—{e-'— - e-'—) = 2(1 - e'^^O 



= 4ismi^-^je^'^. (13) 
o If r(x) = cclw, its contribution to Wo is: 

( 1 ) fvj fwj fzj fwj ~l" fvj fzj fzj fwj + ( 1 ) fzj fvj ( ) /f j fwj — 

~ fwjfzjfvji fwj ~l" fzj ~l" fvj) 

= fwj fz, fv, {Vzj ) (by definition of the vectors /) 
= 2U,fv, = 2e*^(e-^ - e"*^) = 2(1 - e^^^ 
= -4isin(^)e4. (14) 
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Combining this, we deduce the contribution of a decoration x. Note that in computa- 
tions below, we use the fact that ^ = f • 

o If r(x) = cw, the contribution of a decoration x to Wo is: 

n 2i sin(l) I 2'^-i = _22*^-ii<ix+i "Q ^^^(^^'^ . 
o If r(x) = cclw, the contribution of a decoration x to Wo is: 

f](-l)2zsin(|)e4 j (-2)<^-i = -2^'^-- V^+i f] sin(|) . 

We deduce that the contribution of a decoration is in fact independent of its type. Let 
us denote by 

N = V{Gi), M = E{Gi). 

Then, since F is an OCRSF of Gi, it has N edges. Taking the product over ah long 
edges of F yields a contribution: 

(_2)^i^ J] sin^^y, (15) 
(x,y)eF 

where the contribution is in fact independent of the orientation of the edges. Observing 
that X^xey(Gi) ~ 2M, and taking the product over all decorations of Qi yields a 
contribution: 

(_l)iV24M-iV.2M+iV -Q sin2(^). (16) 
xy6£(Gi) 

Taking the product of Equations (15) and (16), gives: 



Wo = 



xyeE(Gi) 



n sinewy. (17) 

xyeF 



Computation of Wfc, fc > 1 . Recall that J^'(^'^*)'{^n'-'^*fc>(gi, g*) consists of all 
pairs {F,F*) of dual essential OCRSFs of Qi and Qi, obtained by performing essential 
moves adding the edge +ejj. or-ej^., on pairs of dual OCRSFs of J"^'*-^'^ '^'''-^\Qi,GI). 
Using Lemma 22, characterizing the set ^-''(^i,^*), and returning to the def- 

inition of essential moves, see Figure 22, we deduce by induction that, for every 
k E {0, - ■ ■ , m} and for every {ii, - ■ ■ ,ik) £ </fe, the restriction to a decoration x of Gi of 
the first component F of a pair {F,F*) of dual OCRSFs oi J^'^^'^'^'^^'i'- '^'k}{g^,g*) 
contains: 



1. all inner edges with the orientation induced by the type r(x), 
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2. any of the three possible 2-edge configurations at triangles of non-root vertices, 

3. the edge {wi,vi) (resp. {zi,vi)) at the triangle ti(x), if r(x) = cw (resp. r(x) = 

cdw), 

4. one or two of the 1-edge configurations at triangles of other root vertices, depend- 
ing on whether one or two essential moves can be performed. 

Let us compute the ratio yp^- To simplify notations, we denote by xy the edge ej^,, 
by Vk{'x.)v£{y) the corresponding long edge of Qi, and by 9 the corresponding rhombus 
half-angle. 

The edge e,^ is absent in OCRSFs of ^i*), so that Vk and are 

non-root vertices. By Point 3., this implies that the triangles and each contain any 
of the three possible 2-edge configurations. Using equations (13) and (14), we deduce 
that the contribution to Wk-i of the absent edge and of the two incident triangles is: 

• If t(x) = cw and r(y) = cw, 

(4zsin(^)e-*f )2 = -16sin2(^)e-^^ 

• If r(x) = cw and r(y) = cclw, or r(x) = cclw and T(y) = cw, 

16sin^(^). 

• If r(x) = cclw and r(y) = cclw. 

We now describe what happens when performing an essential move adding the long 
edge {vk,ve) or {vi,Vk)- Suppose that the dual edge in F* is oriented as in Case A of 
Figure 22. We give all details for the first case only, since others are similar. 

• If r(x) = cw and r(y) = cw. When the essential move adds the edge VkVi in 
either of the two directions, one of Vk or V£ is a non-root vertex, so that the 
contribution of or is: 

4zsm(^-je 2. 

By Equation (12), the contribution of the edge {vk,vi) is independent of its ori- 
entation and is equal to: 

—2i sm.9. 
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Suppose that the essential move adds the long edge {vk,Vi). Then, by Case A 
of Figure 22, there are two possible configurations at the triangle tk, and the 
contribution is: 

fwkfvk^Wk,Vk + fwkfzk^WkjZk — fwkfvk ~ fwkfzk — fwkfwk — 1- 

Suppose that the essential move adds the long edge {vi, Vk)- Then, using symme- 
tries, we are in Case B of Figure 22. Thus there is one possible configuration at 
the triangle t^, and the contribution is: 

As a consequence, the contribution to Wfc of the edge added in one of the two 
possible directions and of the two incidents triangles is 

8sin(^^)e-^isin^(l + e-*^) = 16 sin (^^) e"*'^ sin ^ cos (^^). 
If r(x) = cw and T(y) = cclw. In a similar way, the contribution to Wk is: 
-4isin(|^)e'i^ (-2zsin^)(l) + ^4zsin(^^)e-^i^ (-2zsin0)(-l) 
= — 16sin^-j sin 6* cos ^-j. 
If t(x) = cdw and r(y) = cw, the contribution to Wfc is: 

^4zsin(^)e-^ij (-2isine)(-e*^) + (^-4isin(^)e^i^ (-2isine)(e-^^) 
= — 16sin^-j sin ^ cos ^-j. 
If t(x) = cdw and r(y) = cclw, the contribution to Wk is: 

-4isin(^^)e*ij (-2isin^)(-e'^) + (^Aisui{^e'i^ (-2isine)(-l) 
= lesinf^^e*" sin cos 



Moreover, by Points 1. — 4. we know that the contribution of all other long edges, inner 
edges and triangles is the same in Wk and Wk-i- Thus, from the above computations 
we deduce that, independently of the type of the decorations x and y: 

Wk _ 16sin(|) sin6'cos(|) 



Wfc-i 16sin2(|) 

a 

= -sine cot (-)= -(1 + COS e). (18) 
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Let us now introduce the notation for the rhombus half-angle of the edge ej^. . Then, 
from Formula (18), we deduce that independently of the type of the vertices: 

k 

Wk = mll{-ii + cosei.)). 

Conclusion . Independently of the type of the vertices, we have: 

m m k 

E E >v. = >voE E n(-(i+^°«^^.)) 



k=0 (ii,- ,ik)eJk 



k=0 {ii,- ,ik)&Jk 3='^ 
m 

Wol[[l-{l + cosei)] 

i=l 

m 



\M-N 



JJ cos^xy, (changing notations) 

xy€E(Gi)\F 



Wo(-l 



xyG£;(Gi) 



xy6F 



Jl cos 6'xy 

xyGE(Gi)\F 



2^^ n Sin2(^)cos0: 

xygB(Gi) 



xy 



JJ tan^xy 

xyeF 



(by 



Summing over the 2^ possible types for the vertices of Gi yields Theorem 30. 
Corollary 31 ([BdTlOb, BdTlOa]). 

Pdimer(z,w)=(2l^(«l)l J\ [cot^ (^) - 1]) PLap(z, w). 

xyeBCGi) 

Proof. Prom Theorem 29 and 30, we deduce that: 

^'dimer(z,w) = CPLap(z,w). 

Moreover, by Equation (5), 

^'dimer(z,w)= ( \{ Pdimer(z, w). 

xey(gi) 

By definition of the vector /, for every decoration x and for every k G {I,-- - ,dx}: 

l/«Ife(x) I = l/zfc(x) I = 1 

|/.,(x)P = 2(l-cose) = 4sin2(^), 



□ 
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where 9 is the rhombus half-angle of the long edge incident to t'fc(x). As a consequence: 



n i/.p=2^^ n 

xeViQi) xyeE(Gi) 

We conclude: 

^'dimer(z,w) = (— — ) ^Lap (z, w) 

^llxeViQi) \Jx\ ^ 

2" n ^^U.p(^.w) 

xyGE(Gi) Sm I 



(2^ n [-t^(%)-l])^Lap(z,w). 
xyeE(Gi) 



□ 
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